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In this paper, the holographic dark energy (HDE) model, where the future event horizon is taken 
as an IR cut-off, is confronted by using currently available cosmic observational data sets which 
include type la supernovae, baryon acoustic oscillation and cosmic microwave background radiation 
from full information of WMAP-7yr. Via the Markov Chain Monte Carlo method, we obtain the 
values of model parameter c = 0.696^o'o737-o 132-0 190 with 1,2,3(7 regions. Therefore one can 
conclude that at lest 3<r level the future Universe will be dominated by phantom like dark energy. It 
is not consistent with positive energy condition, however this condition must be satisfied to derive 
the holographic bound. It implies that the current cosmic observational data points disfavor the 
HDE model. 
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I. INTRODUCTION 

In 1998, two supernovae teams discovered that our 
Universe is undergoing an accelerated expansion [l|, 
From that time on, a flood of models have been pre- 
sented to explain the accelerated expansion phase. For 
the reviews about the accelerated expansion and dark 
energy, please see 043 ■ A natural candidate to dark en- 
ergy is the cosmological constant (CC) which was first 
introduced by Einstein to realize a static universe about 
a century ago. In the context of quantum field theory 
(QFT), CC has relations with the vacuum or zero point 
energy density of quantum fields, via 
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where A 3> m is a UV cut-off. To balance an assumed UV 
cut-off A and the observational smallness of CC, tremen- 
dous fine tuning is required. This is the so-called cos- 
mological constant problem Q. QFT in curved space- 
time leads to an infinite effective action or vacuum ex- 
pectation values (VEV) of the energy-momentum ten- 
sors of the fields. A renormalization treatment can yield 
a scale-dependent or running CC and a running New- 
ton constant. Though, the absolute values can not be 
calculated, the change with respect to the renormaliza- 
tion scale can be calculated via renormalization group 
equations (RGEs) originating in QFT [l(| and quantum 
gravity |ll|. In the cosmological context, it is reasonable 
to identify the renormalization scale with some charac- 
teristic scales of cosmology. The related investigation 
can be found in [l2j], where the renormalization scale p 
was given by the Hubble scale H, the inverse radius 
of the cosmological event horizon and the inverse radius 
T _1 of the particle horizon. 



Another interesting viewpoints about the relations be- 
tween vacuum energy density and a cosmological scale 
come from the so-called holographic principle in cosmol- 
ogy. According to the holographic principle, the number 
of degrees of freedom in a bounded system should be fi- 
nite and has relations with the area of its boundary. By 
applying the principle to cosmology, one can obtain the 
upper bound of the entropy contained in the universe. 
For a system with size L and UV cut-off A without de- 
caying into a black hole, it is required that the total en- 
ergy in a region of size L should not exceed the mass of 
a black hole of the same size, thus L 3 p\ < LM 2 ^ The 
largest L allowed is the one saturating this inequality, 
thus p\ — 3c 2 M 2 t L^ 2 , where c is a numerical constant 
and M p i is the reduced Planck Mass M~ 2 = 8nG. It 
just means a duality between UV cut-off and IR cut-off. 
The UV cut-off is related to the vacuum energy, and IR 
cut-off is related to the large scale of the universe, for 
example Hubble horizon, event horizon or particle hori- 
zon as discussed by [UflBl- Then in this context the 
so-called fine tuning problem is removed naturally. 

In the paper , the author proposed a dark energy 
model, dubbed holographic dark energy (HDE), where 
the future event horizon 
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was taken as the IR cut-off L. This horizon is the bound- 
ary of the volume a fixed observer may eventually ob- 
serve. One is to formulate a theory regarding a fixed 
observer within this horizon. As pointed out in (To| . it 
can reveal the dynamic nature of the vacuum energy and 
provide a solution to the fine tuning and cosmic coinci- 
dence problem. 

The concrete value of parameter c is very important 
to determine the properties of HDE. For example in the 
epoch when the HDE becomes to dominate the energy 
density, i.e. fth ~ 1, if c > 1, c = 1 and c < 1 are 
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respected, the HDE behaves like quintessence, cosmolog- 
ical constant and phantom respectively. One can clearly 
see that in the following Eq. (fT0|) by setting fi/, = 1. 
In the past years, the HDE model was confronted with 
cosmic observations. In Ref. [l6j], the authors used the 
angular scale of the acoustic oscillation from the BOOM- 
RANGE and WMAP data on CMB to constrain HDE 
model. Also the authors of Ref. [l7| used the CMB 
shift parameter R to constrain HDE model. For the re- 
cent result please see fl8l | where the best fit value of pa- 
rameter c = 0.65^o'o8 m l* 7 region was given. In that 
paper, about the CMB information only the CMB shift 
parameter R, acoustic scale Ia and the redshift z* at 
the decoupling epoch of photons obtained from WMAP- 
7yr in ACDM model were employed. The main problem 
comes from the fact that the values of the CMB shift 
parameters used in the literature for example [l8[, de- 
pend on the cosmological model, here the ACDM model. 
So, when one uses these derived data points to constrain 
other cosmological models, the so-called circular prob- 
lem, i.e. dependence on cosmological models, would be 
committed. And the circular problem is what we try to 
avoid in the cosmological constraint issue. Therefore one 
should constrain the cosmological model in a consistent 
way by using cosmic observational data points which do 
not depend on any cosmological model. However to use 
full CMB data sets from WMAP7-yr to constrain the 
HDE model is still missing in the literature. In this pa- 
per to fill up the gap, we will use the full information from 
WMAP-7yr data sets in stead of CMB shift parameters. 
Furthermore, the CMB shift parameters contain infor- 
mation around the last scattering surface of background 
photons where the contribution from dark energy may 
not be neglected as CC case due to a nontrivial equation 
of state (EoSjof dark energy, for example early dark en- 
ergy model In this situation, it is dangerous to use 
theses derived data points from ACDM model to con- 
strain other cosmological model due to much departure 
from ACDM model. At late epoch when the dark en- 
ergy dominates, the gravitational potentials decay and 
the late integrated Sachs- Wolfe (ISW) arises. It affects 
to the anisotropy power spectra of CMB at large scale 
(low I parts). The ISW effect is sensitive to the proper- 
ties of dark energy. Apparently, the late ISW effect is not 
expressed by the CMB shift parameters. So one can ex- 
pect a tight constraint to cosmological model parameter 
space when the late ISW effect is included. 

Base on these points mentioned above, in this paper we 
will not use the derived CMB shift parameters but the full 
information from WMAP-7yr data points to constrain 
the HDE model. The important is that there doesn't 
exist any circular problem. Due to the full information 
from WMAP-7yr is used, one can expect to obtain a tight 
constraint. At low redshift regions, the type Ia supernova 
(SN Ia) and baryon acoustic oscillations (BAO) are also 
included as geometric measures. 

This paper is structured as follows. In section UH we 
give a brief review of the HDE model where the radiation 



is included and the future event horizon adopted as an IR 
cut-off. The scalar perturbation evolution equations for 
a spatially flat FRW Universe will also be presented. In 
section IIII1 the constraint methodology and results will 
be presented. We give a summary in section HVl 



II. BACKGROUND AND PERTURBATION 
EVOLUTION EQUATIONS 

Now, as is done in fl5j , the event horizon R e h taken as 
the IR cut-off. Then, the HDE is written as 
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R lh 



(3) 



and the Friedmann equation in a spatially flat FRW Uni- 
verse can be written as 

H 2 = H 2 (O-ofl- 4 + n b0 a- 3 + n c0 a- 3 ) + n h H 2 , (4) 

where f2i = nj. fi H g are dimensionless energy densities 

pi 

for radiation, baryon, cold dark matter and HDE respec- 
tively. For convenience, the scale factor a has been nor- 
malized to do = 1. Combining Eq. ([3]) and Eq. ©, one 
has 
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where the definition of dimensionless energy density pa- 
rameter for HDE is used. From Eq. (J4j , one obtains 
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Inserting the above equation into Eq. ([5]), one has 
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where E(x) — \/fl r ae~ 2x + flboe~ x + fl c oe~ x and x = 
In a. Taking derivative with respect to x = In a from 
both sides of the above equation, one has the differential 
equation of 0,^ 



n' h = -2n h (i - n h ) 
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where ' denotes the derivative with respect to x = In a. 
This equation with initial condition ttho = 1 — f2 r o — 
flbo — r2 c o describes the evolution of dimensionless energy 
density of dark energy. Then the background evolution 
will be obtained by solving this differential equation and 
Friedmann equation. When the baryon and radiation 
components are omitted, one has E'(x)/E(x) = —1/2. 
Then the evolution equation for O/, [1 51 ] 



n' h = n h {i-n h ) + 
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is recovered. From the conservation equation of the HDE 
Ph + 3H(p h +p h ) = 0, one has the equation of state (EoS) 
of holographic dark energy 
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where w h = p h /ph- 

In this paper, we take the HDE as a perfect fluid 
having EoS (fTU)) . Considering the perturbation in 
the synchronous gauge and the conservation of energy- 
momentum tensor = 0, one has the perturbation 
equations of density contrast and velocity divergence for 
the HDE 



h = -(1 + w h )(6 d + -)- m{c 2 s - w h )S h (11) 
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following the notations of Ma and Bertschinger [20| . For 



the gauge read y fo rmalism about the perturbation the- 
ory, please see For the HDE, we assume the shear 
perturbation ah = 0. In our calculations, the adiabatic 
initial conditions will be taken. 



III. METHODOLOGY AND CONSTRAINT 
RESULTS 

Before constraining the model parameter space, we 
study the effects of model parameter c to the CMB power 
spectrum by setting different values of c and fixing the 
other relevant cosmological mod el p arameters. At fist, 
we modify the CAMB package [22|], which is the code 
for calculating the CMB power spectrum and is publicly 
available, to include the HDE. The differential equation 
([8|) was solved by taking QhO = 1 — ^ m o — ^Vo a s initial 
condition. We borrow and fix the other relevant cosmo- 
logical parameters values from WMAP7, the resultant 
CMB power spectra are shown in Fig. [TJ for different 
values of c, where h = 0.72, uj c = 0.112, uj b = 0.0226, 
n., = 0.96 are fixed. 
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From Fig. [TJ one can see that the main contribution 
to CMB power spectra from the HDE comes from the 
evolution of EoS and its domination at the late epoch 
through late ISW effect. At early epoch the contribution 
to CMB power spectra is small due to small energy den- 
sity ratio of the HDE, that can be understood from the 
EoS of the HDE (fTU]) which seems close to —1/3 before 
the last scattering surface. It means that a tight con- 
straint would be obtained when CMB information from 
large scale is included. But unfortunately, the CMB data 
points are dominated by the so-called cosmic variance at 
large scale, i.e. low I part. 

We perform a global fitting to the model parameter 
space by using the Markov Chain Monte Carlo (MCMC) 
method. We modified publicly available cosmoMC 
package Q to include the HDE in the CAMB [22| code 
which is used to calculate the theoretical CMB power 
spectrum. The following 7-dimensional parameter space 
is adopted 

P = {cu b , u c , G s , r, c, n s , log[10 10 A s ]} (13) 

where oj b = il b h 2 and oj c = Q c h 2 are the physical den- 
sity of baryon and cold dark matter respectively, &s 
(multiplied by 100) is the ration of the sound horizon 
and angular diameter distance, r is the optical depth, c 
is the newly added model parameter related to HDE, 
n s is scalar spectral index, A s is the amplitude of of 
the initial power spectrum. The pivot scale of the ini- 
tial scalar power spectrum fc s o = 0.05Mpc _1 is used 
in this paper. The following priors to model param- 



eters are adopted: uj b G [0.005,0.1], cj c G [0.01,0.99], 
6 S G [0.5,10], t G [0.01,0.8], c G [0,1], n s G [0.5,1.5], 
log[10 10 ^4 s ] G [2.7,4]. Furthermore, the hard coded prior 
on the comic age lOGyr < to < 20Gyr is also imposed. 
Also, the physical baryon density u b = 0.022 ±0.002 0] 
from big bang nucleosynthesis and new Hubble constant 
H = 74.2 ± S.ekms^Mpc" 1 [H are adopted. 

To get the distribution of parameters, we calculate the 
total likelihood C cx e~ x / 2 , where \ 2 is given as 

X 2 = XCMB + X 2 BA0+XSN- ( 14 ) 



The 557 Union2 data ^26| with systematic errors and 
BAO [i2l are used to constrain the background evolu- 
tion, for the detailed description please see Refs. [29j . 
SN la is used as standard candle. And BAO is used as 
standard ruler. At the last scattering of CMB radiation, 
the acoustic oscillation in the baryon-photon fluid was 
frozen and imprinted their signature on the matter distri- 
bution. The characterized scale of BAO in the observed 
galaxy power spectrum is determined by the comoving 
sound horizon at drag epoch Zd which is shortly after 
photon decoupling and defined as the redshift at which 
the drag optical depth Td equals one. To calculate r s (zd), 
one needs to know the redshift za at the drag epoch and 
its corresponding sound horizon. The baryon drag epoch 
redshift Zd can be obtained from its definition and be 
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FIG. 1. The CMB Cj T power spectrum v.s. multiple moment I and the corresponding ratio to ACDM model for different 
values of model parameter c, where h — 0.72, uj c — 0.II2, uji — 0.0226, n 3 = 0.96 are fixed. The main difference appears at low 
(I < 20) multipole momentum parts which correspond to large scale. 



calculated numerically [3C 
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where R = 3pfc/4p 7 , <jt is the Thomson cross-section and 
x e (z) is the fraction of free electrons. Then the sound 
horizon is 



rvi z d) 

r s (z d ) = / d-qc s (l + z). 
Jo 



(16) 



where c s = l/y/3(l + R) is the sound speed. Also, to 
obtain unbiased parameter and error estimates, we use 
the substitution 130] 



, r s (z d ) 

d z —, — rr s (z d ), 
r s (z d ) 



(17) 



where d z = r s (z d )/ 'Dv(z), f s is evaluated for the fiducial 
cosmology of Ref. [23|, and z d is redshift of drag epoch 
obtained by using the fitting formula [3l| for the fidu- 
cial cosmology. Here D v (z) = [(1 + z) 2 D\cz/ 'H(z)} 1 ^ 3 
is the 'volume distance' with the angular diameter dis- 
tance Da- In this paper, for BAO information, the SDSS 
data points from [271 ] and the WiggleZ data points [28[ 
are used. For CMB data set, the temperature power 
spectrum from WMAP 7-year data [32| are employed as 
dynamic constraint. 

After running 8 independent chains and checking the 
convergence to stop sampling when the worst e-values 
[the variance(mean)/mean(variance) of 1/2 chains] R—l 
is of the order 0.01, the global fitting results are summa- 
rized in Table fl] and Figure [2j We find that the values 
of c in 1,2,3c regions are c = 0.696^ 



H0. 0736+0. 159+0.264 
'-0. 0737-0. 132-0,. 



which are compatible with the result obtained in 
And the error bars are slightly smaller than that in 
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Prameters Mean with errors Best fit 



Q,DMh 

9 

T 

C 

n s 

logflO 10 ^ 

V h 
Age/Gyr 

H 



n n99R+°' 000536 + ' 00108 

U.UZZO_p .000546-0 .00107 
n i i n +0. 00422+0. 00866 
U.llU_o. Q0414-0. 00838 
h0. 00268+0. 00513 
-0.00263-0.00520 

n n«Q9+°' 00635 + ' 0251 

u - u ° yz -0. 00719-0. 0226 
fl CQR+0. 0736+0. 159 
U.O»O_ .0737-0.132 
fl Q71 +0.0135+0.0262 
U.y / 1_ . 0134-0. 0255 
o nS no+0-0340+0.0704 
O.UOUZ_ .0337-0 .0646 
n 7 9 7 +0.0172+0.0323 
u - lz ' -0.0171-0.0349 
1 878+°' HO+0. 221 
lo -°'°-0. 113-0. 216 
fl 970+0. 0171+0. 0349 
u - z '°-0. 0172-0. 0323 
10 ROQ+1.189+2.352 
±u - ozo -l. 161-2. 229 
RQ Q1 e+1.867+3.635 
Qg-gJ-D-1. 844-3. 585 



0.0226 
0.112 

1.0400 
0.0865 

0.643 

0.971 
3.0805 

0.733 
13.839 

0.267 
10.468 
71.0657 



TABLE I. The mean values of model parameters with la and 
2a errors, where WMAP 7-year, SN Union2 and BAO data 
sets are used. 



This result shows that the values of c are less than 1 at 
3<t level. Correspondingly, we plot the CMB power spec- 
trum for the mean values estimated from MCMC anal- 
ysis in Figure [31 where the observational data points of 
7-year WMAP with uncertainties are also included. As 
comparisons, the CMB temperature power spectrum for 
ACDM with same cosmic observational data sets combi- 
nation and results from [32| are also shown. One can see 
that the CMB power spectrum for the HDE is well inside 
the error bars of the binned measurements from WMAP 
7-year results and matches to ACDM model very well. 
This strongly implies that current cosmic observational 
data combinations of SN Union2, BAO and WMAP 7- 
year can not almost discriminate the HDE from ACDM 
model. In this sense, the HDE is a competitive model of 
dark energy. 




FIG. 2. The ID marginalized distribution on individual parameters and 2D contours with 68% C.L. and 95% C.L. by using 
CMB+BAO+SN data points. The shade regions show the mean likelihood of the samples. 
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FIG. 3. The CMB C{ 1 power spectrum v.s. multiple moment I, where the grey dots with error bars denote the observed data 
with their corresponding uncertainties from WMAP 7-year results, the green solid line is for HDE with mean values as shown 
in Table U the red dotted line is for ACDM model with mean values taken from [32| with WMAP+BAO+ifo constraint results. 



We can show the evolutions of dimensionless energy 
density of the HDE and the EoS of the HDE with respect 
to redshift z in Fig. 21 one can solve the differential 
equation ([7]) by setting the initial conditions f2/jo and 



model parameters f! m o = flbo + Q c o and c to their central 
values as listed in Table U here fi r o is derived from the 
temperature of CMB background T = 2.726X. 
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o.o — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 -1.1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 

2 4 6 8 10 02468 10 

z z 

FIG. 4. The evolutions of dimensionless energy density of the HDE (left panel) and EoS of the HDE (right panel) with respect 
to redshift z, where the central values of model parameters are taken as listed in Table [I] 



From the Fig. |4j one can read off that the HDE is phan- 
tom like at late time and quintessence like at early epoch. 



IV. SUMMARY AND DISCUSSION 

In this paper, we perform a global fitting to the HDE 
model by using MCMC method with the combination of 
the full information from CMB, BAO and SN la data 
points. We analyze the effect of model parameter c to 
the CMB power spectra by taking its different values. We 
find out that c mainly contributes to CMB power spectra 
at large scale where the so-called cosmic variance is dom- 
inated. However, we still obtain a relative tight constrain 
to model parameter space where the values of c with 
1, 2, 3a regions are achieved: c = 0.696±g;^lf±g:i||±g;?|^. 
So at 3cr level, c is a number less than 1. As a comparison 
to the previous results obtained in [ijj], a slightly tight 
constraint is achieved in this paper. More important is 
that a consistent constraint is carried out instead of us- 
ing the derived CMB shift parameters. Based on this 
result, one can conclude that at least 3cr level the future 
Universe will be dominated by phantom like dark energy. 
Then the positive energy condition is not respected, how- 
ever this condition must be satisfied to derive the holo- 
graphic bound. So the current cosmic observational data 
points disfavor the HDE model based on our analysis. 

Based on the holographic principle, other HDE models 
were proposed in the literature where different IR cut- 



offs were adopted. For example, the Hubble horizon as 
IR cut-off was discussed when the model parameter c is 
constant [3, 0j| or time variable [33[ . In Ref. [1H , the 
authors took the Ricci scalar as the IR cut-off and named 
it Ricci dark energy, the corresponding cosmic observa- 
tional constraint to this model could be found in Ref. 
[35l ]. Interestingly, Cai, et. al. found out that the holo- 
graphic Ricci dark energy had relations with the causal 
connection scale R7?n = Max(H + 2H 2 , —H) for a spa- 
tially flat Univers e 13611 . Based on the idea that gravity as 
an entropic force [37J , a similar DE density was given in 
[38| where a linear combination of H 2 and H was also pre- 
sented, see also [3^, H3| ■ And a generalized HDE model 
p R = 3c 2 M 2 t Rf (H 2 / R) and Ph = 3c 2 M 2 l H 2 g(R/ H 2 ) 
were presented in Ref. 41]. To the best of our knowl- 
edge, the full information from CMB has not been used 
to test these models yet. One can apply the method and 
data sets presented in this paper to these models directly. 
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